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Abstract. We analyze the continuous time evolution of a d-dimensional system of N 
self propelled particles subject to a feedback rule inspired by the original Vicsek's one 
[VCB -TTJS] . Interactions among particles are specified by a pairwise potential in such a 
way that the velocity of any given particle is updated to the weighted average velocity 
of all those particles interacting with it, which makes the system non-Hamiltonian. 
The weights are given in terms of the interaction rate function. When the size of the 
system is fixed, we show the existence of an invariant manifold in the phase space 
and prove its exponential asymptotic stability. In the kinetic limit we show that the 
particle density satisfies a Boltzmann-Vlasov equation under suitable conditions on the 
interaction. We study the qualitative behaviour of the solution and we show that the 
Boltzmann-Vlasov entropy is strictly decreasing. 



1. Introduction 

The analysis of the evolution of a network of a large number of coordinated self 
propelled particles (agents) is a sub discipline of control theory which is has seen a 
rapid development during the last decade jBDTl Ml iJLMl ICKFLl IB-NVRl [CSllCHDBj . 
This is due to its several potential application in understanding the emerging of collective 
behavior in biological systems (for example fish schools and bird flocks), computer science 
jEtlBPT] . engineering jJLMl ICSl ICiHDB] . economy [DY] and social sciences [Wl ICKFLl 
IB-NVRj . To model the particle self-organized behavior one assigns to any particle a 
simple communication/interaction rule in order for the whole system to dynamically 
reproduce, in a given regime of the model's parameters, specific phase space patterns. 

The emergence of phase space patterns persistent in time described by a large con- 
nected cluster of coherently moving particles is called flocking or swarming (also school- 
ing or herd behavior). Basic models of flocking behavior generally follow three simple 
rules: 1) separation, that is to avoid crowding neighbors (usually modeled by short range 
repulsion interactions); 2) alignment, i.e. to steer towards average heading of neighbors; 
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3) cohesion, i.e. to steer towards the average position of neighbors (usually modeled by 
long range attraction interactions). 

The seminal work in the direction of modeling flocking behavior is the one of Vicsek 
et al. |VCB- JCS] . They proposed a model of N interacting particles located on a 
2-dimensional torus of diameter D. The velocity of each given particle belongs to the 
unit circle and at each time step its direction is updated at the empirical average of 
the velocity's directions of all the particles lying in a neighborhood of radius 1 from 
the given one, including itself, plus a random perturbation. Particles positions are 
then updated according to their velocity. Computer simulations proved that, when the 
particle density is sufficiently high and the noise intensity sufficiently small, the 
distribution of the velocities of the particles concentrates around the velocity of the 
barycenter of the system, although this is not a quantity preserved by the dynamics. 

We propose a simple model of continuous time noiseless multi-agent evolution closely 
inspired to the original Vicsek's one. The particles interact (communicate) with each 
other trough a pairwise interaction function, which can be chosen to have the shape of 
an ordinary interaction potential, in such a way that the velocity of any given particle 
is updated to the weighted average velocity of all those particles communicating with 
it, with a weight given in terms of the communication rate function. This makes the 
system non-Hamiltonian. For what concerns flocking behaviour our model takes into 
account allignment and cohesion, but violates the separation rule since the particles can 
overlap. 

We prove for such model two type of results. First, we analyze the N particle dynamics 
in M d . We show that there exists an invariant manifold in the phase space and prove 
exponential asymptotic stability of the invariant manifold when the initial conditions for 
particles dynamics are suitably chosen. This implies that the system, under the chosen 
initial conditions, will reach a state of flocking. Then, we will study the mean-field limit 
(N — > oo) of the system and prove that the particle density satisfies a Boltzman-Vlasov 
equation when the particles are confined on a torus and subject to a short-range potential 
of Gaussian type. The same result holds in M. d when the interaction among particles is 
given by a suitable regularization of a finite range potential. We further show that the 
Boltzmann-Vlasov entropy is strictly decreasing. As a consequence, one can argue that, 
even if the initial distribution of the particles is absolutely continuous w.r.t. Lebesgue 
measure, the limit density distribution is singular w.r.t. Lebesgue measure. 

A continuous time version of Vicsek's model, as well as its stochastic counterpart 
driven by the Brownian motion, has been proposed in |DM] and the corresponding kinetic 
equations heuristically derived and studied. In fact, at present time, to our knowledge, 
a rigorous derivation and analysis of Vicsek's model kinetics, as well as hydrodynamics, 
is lacking. 

Another basic model for flocking is the Cucker-Smale one |CSj . In this and related 
models [DQBCj IAIR] the variation in time of the momentum of a given particle is the 
weighted sum of the differences between the particle's momentum and those of the other 
system's components, with weights depending of the relative distances among particles 
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divided by the total number of particles N. It is worth notice that, for all these models, 
the interaction among two given particles is of order 1/N, therefore when the size of 
the system becomes large, particles tend to decorrelate. On the contrary, in the original 
Vicsek's model, the interaction between a given couple of particles is of order one. 
Moreover, Cucker-Smale dynamics preserves the velocity of the barycenter, which is not 
the case for Vicsek's. As a matter of fact, we believe that the order of the interaction 
with respect to the size of the system is the peculiar feature distinguishing Vicsek's 
from Cucker-Smale algorithm. Therefore, in our opinion, variants of the Cucker-Smale 
momenta updating rule taking into account only the differences among the directions 
of the momenta of the particles, rather than those of the momenta as vectors, are 
somewhat improperly ascribed to variants of the Vicsek's model |BCC2j . Cuker-Smale 
and related models have been more deeply investigated in the mathematical literature 
and their mean-field limit equations rigorously derived and studied in |HTl IHL| ICFRT} 
ICCR] lAIRj in the noiseless case and in [BCClj IBCC2] in the stochastic case driven by 
Brownian motion. Moreover, the hydrodynamics equations for these models have also 
been rigorously studied but formally derived [HTt ICDP[ ICCR] . 

The plan of the paper is the following. In Section 2 we describe the model and set the 
notations. In Section 3 we analyze the system when the number of particles is fixed. In 
Section 4 we analyze the system when the number of particles goes to infinity. In the 
appendix we collect proofs of results used along the previous sections. 

Acknowledgements. Enza Orlandi thanks Carlangelo Liverani for useful discussions. 
Michele Gianfelice thanks Fabio Fagnani and Marco Isopi for interesting discussions on 
the subject. 



2. Description of the model and notation 

2.1. Notations. Given x G M. d ,d > 1, we denote by x l its i-th component, % — 1, ..,d, 
with respect to the canonical basis (ei, .., ed) ■ For any x, y G M. d we set x ■ y := Yut=i x% V l 
to be the scalar product between x and y. Hence, we denote by \x\ := \Jx ■ x the 
associated Euclidean norm and by B r (x) := {y G M d : \y — x\ < r} the ball of radius 
r > centered in x and B r := B r (0) . Furthermore we set H^H^ := maxj = i r . i( j \x % \ . 

Given an integer N > 2, let x := (xi, .., Xjy) G WL Nd . We denote by x • y := Yji=i x,i ■ yi 
and |x| := v /x T ~y respectively the scalar product and the Euclidean norm in M. Nd . We 
also denote by B r (x) := {y G M> Nd : |y — x| < r} the ball of radius r > centered in x. 

Partial derivative w.r.t. any component x l of x G M. d will be denoted by d x i, so that 
V x stands for (d x i, .., d x a) while, for any q G M. Nd , we set V q := (V 91 , .., V gjv ) . 

Moreover, we denote by C n (M) the Banach space of linear operators from M. n to itself 
and by ||-|| and IHI^ the operator norm induced by respectively the Euclidean and the 
supremum norm. In particular I n , n G C n (R) denote respectively the identity and the 
null operator. 
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2.2. The model. Let N > 2 be an integer. We consider N particles of unitary mass in 
M. d evolving according the equations: 



where, for i = 1, ..,N, (qi,Pi) G R d x R d , (c?°,p°) are the initial conditions and £/#(•) is 
the two particles interaction. We assume that Ur(-) is a sufficiently smooth, spherically 
symmetric positive function, actually a probability density concentrating the mass in 
the ball of radius R centered at zero, such that sup^^d Ur(x) = Ur(0). This implies 
that the denominator in the second equation of (12. ip is always strictly positive. Other 
assumptions on Ur will be done in Section 3, where the kinetic limit as iV — > oo will 
be considered. The vector field in (12.11) is Lipschitz, therefore the existence and the 
uniqueness of the solution is granted at least for short time. Since the vector field 
increases at most linearly in (q, p) the solution exists and it is unique for all t > 0. 

2.2.1. Communication graph. We recall some basic definition of graph theory and its 
applications to Markov chains used in the next section. We refer the reader to basic 
textbooks such as [B] and [St] for an account on this subject. 

A directed graph Q is a ordered pair of sets (V, £) where V is a finite set called set 
of vertices and £ C V x V is called set of edges or bonds. Q 1 = (V',£') such that 
V C V and £' C (V' x V') fl £ is said to be a subgraph of Q and this property is 
denoted by Q' C Q. Two vertices are said to be adjacents if belong to the same bond, 
so that u, v G V are adjacents w.r.t. b G £ if either b = (u,v) or b = (v,u) . For any 
b G £, let then V& '■= {{u,v} cV:u and v are adjacent w.r.t. b} . Then, if b = (u,v) , 
b is said outgoing from u and ingoing in Let := {6 G £ : 6 = (f , u) , ii G V} be 
the set of bonds outgoing from v. We denote by A^ + (v) := (u 6g£: +V&) C V the closed 
outgoing neighborhood of t> and set, for any W C V, A/" + (W) := U t , e wA/' + (f ) to be 
the closed outgoing neighborhood of W. Given v G V, we set A/"^ (u) := Af + (v) and, 
for k > 2, A/"^ (f ) := A" + (A/^]^ (f )) to be the outgoing k -neighborhood of v. Given 
two vertices u and v , u is said to communicate with t> if there exists > 1 such that 
w G A/^ (f ) . Therefore, u, v G V are said to be connected if one communicates with the 
other. In fact, since if w G A/^ (f ) for some k > 1, then u G A/" ; + (f ) , Wl > k, for w and v 
to be connected there must be hi, ki > 1 such that w G A/S (v) and u G A/j^ (u) , that is 
m G A/"^ vA , 2 («),»£ ■^k 1 wk 2 ( M ) • ^ ^ s then said to be strongly connected if any two distinct 
vertices are connected. The maximal connected subgraphs of Q are called components 
oig. 

An example of directed graph is the one which can be associated to a Markov chain. 
In this case, V coincides with the set of states of the chain and, denoting by P the 
transition matrix associated to the chain £ = £ (P) := {(u, v) G V x V : P UjV > 0} . 
Then, the directed graph associated to the Markov chain with transition matrix P is 




(2.1) 
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denoted by Q (P) . Hence, the chain and therefore P are said to be irreducible if and 
only if Q (P) is strongly connected. 

Definition 2.1. Given any particles configuration q G R Nd , we define the geomet- 
ric communication graph of the particle system to be the directed graph <3r (q) : = 
(V (q) , (Br (q)) , where V (q) is the finite subset ofR d associated to q and 

£r (q) := {(<?, q) G V (q) x V (q) :U R (q- q') > 0} . (2.2) 

When considering the particles configuration at a given time t, q (t) , we set Vt ■— 
V(q(f))and<S fl (*):=<Mq(*))- 

Definition 2.2. The system is said to have reached a state of flocking if there exists 
v G R d such that, for any e > 0, 3T e > : Vt > T e , 

• Pi {t) eB e {v)yi = l,..,N; 

• the geometric communication graph (3r (t) is connected and \Vt\ > 2. 

Definition 2.3. The state of the system (q, p) such that qi = ■ ■ ■ = qNiPi = ■ ■ ■ = p^ 
is called rendez-vous state. 

From (12. ip it follows that rendez-vous states belong to the collection of the invariant 
states for the dynamics. 

3. Particle dynamics 

In the following we analyze the evolution of N particles according equations (12. ip . In 
this section N is kept fixed, so we omit in the notation to write explicitly the dependence 
on N. We show that there exists an invariant (N + l)d manifold for evolution (12.11) and 
prove first its stability, see Corollary 13. 5\ and then its asymptotic stability, see Theorem 

Em 

3.1. Stability. We first notice that if the velocities of the particles at time zero are 
bounded, that is, for all i = 1, .., N, p® G B T for some r > 0, then they will lie in B r for 
later times. In fact we have the following result: 

Lemma 3.1. For any i = 1,..,N, assume that Pi{0) G B r . Then, pi(t) G B r , for all 
t > 0. 

Proof. Assume, without loss of generality that r = 1 and that there is a t* such that 
there is at least one Pi(t*) such that = 1 and |pj(£*)| < 1 for j ^ i. Then 

id. Ef=i ggfag) - %•(**)) wr) - p»g*)] • Pig) ^ n 

-~r \Pi(t = — it? < . 3.1 

□ 



Remark 3.2. The result of Lemma \3.1\ holds for any positive smooth interaction Ur, 
regardless of its support. In particular, it holds if Ur does not have compact support. 
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Remark 3.3. The only critical point of the system A2.1\) is (0, 0) . Moreover, if the 
particles at initial time have all the same velocity, that is, for j = 1, ..,N, p°- = v G R d , 
their velocity will remain constant during the evolution and the system describes the 
motion of N non-interacting particles. 

From Remark 13.31 one deduces that the (N + l)d linear manifold 

X = U {vmd} l{v) (3.2) 

where 

X{v) = {(q, p) G R dN x R dN : Pi = v, % = 1, .., N} (3.3) 

is invariant for the evolution (|2.ip . Namely, if the initial data belong to X (v) the particles 
evolve independently one from the other with constant velocity v. 

The next result shows that if at time t = the velocity of the N particles is close to 
its mean velocity vector, then, at any further time t, it will always remain close to the 
mean initial velocity vector. Let Q G C^d be the operator such that 

R Nd 3 x ^Ox G R Nd , (3.4) 

where fix is the vector in R Nd whose component are the vectors (£lx) i = YljLi x j ^ ^ d 
for i = 1,..,N. Notice that by definition Q is the orthogonal projector on 
{x G R Nd : aci = ■ •• = x N ] . 

Theorem 3.4. Let w(t,w°) = (q(t),p(t)) be the solution of \2. II) at time t starting from 



w 



- 'q°, p°) G R 2Nd . Given e > 0, assume that |p° - fip°| < e. T/jen 



|p(*) - Op°| < e, Vt>0. (3.5) 

Proof. We proceed as in the proof of Lemma |3~T1 Let us denote by V{ (t) = Pi (t) — (fip ) i G 
= 1, .., N, and assume that there is a £* so that there is at least one components 
\vi(t*) \ = e and \vj(t*)\ < e for j ^ i. We have 

^W*)| 2 = «i(*)'^i(*) (3-6) 

Ejli - <&(«*)){«i(f ) - )} ■ ^(r) 

< o . 



□ 



Note that, for any u> G 



D2Vci 



dist (w,X) = inf \w — w I = inf |p — p°| = |p — Opl , (3.7) 

u.°6X {p°eR Ard :t« = (q ,p )eX} 

where f2 is the operator defined in (13.41) . From Theorem 13.41 one deduces that the 
invariant manifold X is stable for the evolution (12.11) . 
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Corollary 3.5. For any e > let B(e,X) = \w G M? Nd : dist (w,X) < e} be a neigh- 
borhood of radius e of X. Let w(t,w°) be the solution of (OOP at time t starting from 



w 



q°,p°) G B(e,X). Then 



dist (w(t,w°),X) < 2e, Wt > . 

Proof. By (13.71) we have 

dist (w(t,w°),X) = \p(t) - fip(t)| < \p(t) - ftp ! + \Mp(t) - ^P°| 
By definition of Q, see (13.41) . 

\Qp(t) - fip°| = \n{p{t) - Qp°)\ < \p{t) - Qp°\ . 
Hence, by Theorem I3.4[ 

dist (w(t,w°),X) < 2\p(t) - fip°| < 2e , Vt > . 



(3.8) 

(3.9) 

(3.10) 

(3.11) 
□ 



3.2. Asymptotic stability. Next we show a stronger result. Choosing suitably the ini- 
tial datum, the solution of (12. ip converges exponentially towards the invariant manifold. 
We show that the N particles will not split into non interacting groups and the velocity 
of each particle converges exponentially fast to a velocity vector which is the same for 
all the N particles. In other words, the system will reach a state of flocking as given in 
Definition O 

We rewrite the non linear system (12. ip as follows: 



dt 

q(0) = q°,p(0) = p° 



C(q(f)) 



q(*) 

p(*) 



where 



ONd ^-Nd 

Nd L (q) 



3q^C(q) := 

L(q) :=A(q)-I Nd 
and A(q) is the linear operator valued function so defined 

ai,i(q)Id Oi j2 (q)Id 
a2,i(q)Id fl2,2 (q)Id 
ajv,i(q)Id 



G C 



2Nd 



3 q 



A(q) := 
^fl(g» ~ Qj) 



a liN (q)l d 
a 2 ,7v(q)Id 
dN,N-i(q)h a N:N (q)I d 



1,..,N, 



Remark 3.6. Notice that for q G 



»Nd 



"hi 



[q) = Oij(q + Ox), Vx G 



l,..,iV. 



l,..,iV, 



l,..,iV 



(3.12) 

(3.13) 
(3.14) 

G (K) 
(3.15) 
(3.16) 

(3.17) 
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and 

N 

5>,,-(q) = l. (3.18) 
These two properties are important when studying the spectrum of C (q) for a fixed value 

o/q- 

3.2.1. Spectral Analysis of C (q). Let q e R be fixed. The eigenvalues of C (q) are 
the roots of the characteristic equation 

Det [C (q) - XhNd] = (-\) Nd Det [L (q) - XI Nd ] = 0. (3.19) 

We need then to study the spectrum of L (q) and therefore, by ( 13. 14ft the spectrum of 
A(q). To do this it is convenient to introduce the tensor space R N £g> R d . We denote by 
J 7 the isomorphism 

N d 

R Nd 3 x — >■ J 7 (x) := ^ ^ afe g> e,- eR N ®R d , (3.20) 

i=l i=i 

such that J-"(x)jj = a^, % — 1, .., iV and j = 1, .., d. 

To ease the notation we omit in the following to write the dependence on q if no confu- 
sion arises. We therefore set A := A (q) . One obtains immediately that A : ~R Nd — > R Nd 
acts on ~§i N eg) M d as follows 

A ® I d : R N ® R d — ^M 7V ®M a! , (3.21) 
where, by (13.161) . setting a i; j := ajj(q), 



.4 

Namely, one has that 



a l,l a l,2 ■ ■ ■ a l,N 

a 2,l a 2,2 ■ ■ ■ a 2,N 

a^,i ■ ■ ■ 0>N,N-1 o-n,n 



(3.22) 



(A <8) I d ) 7"(x) = 7 (Ax) . (3.23) 

Furthermore, denoting by C C the spectrum of A, 

E(A) = E(i®I d ) = £(i)£(I d f|. (3.24) 

Since the only eigenvalue of ILj is 1 with multiplicity d, the problem is reduced to study 
the spectrum of A. The matrix A is a (right) stochastic matrix, that is it has non- 
negative entries and, by (13.181) . X^=i a i,j = 1, Vi = 1, .., iV. Then, if it is irreducible one 
can apply the Perron-Frobenius Theorem. We recall that a matrix D 6 £ n (R) with 
non-negative entries is said to be irreducible if there exists an integer m so that _D m has 
strictly positive entries. 



1 If Z :— {z\, .., z n } and W :— {w\, .., w m } are two discrete subsets of C we denote by 

ZW :— {ziWj G C : i = 1, .., n ; j = 1, .., m} . 
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The existence of m G N such that A m (q) has strictly positive entries is related to 
have the graph Q \ A (q)^ associated to A(q) strongly connected. In fact, Q(A(q)) 

is strongly connected if and only if A (q) is irreducible. Clearly this property is also 
shared by &r (q) which is nothing else but a geometric representation of Q{A (q)) in Mr. 
Therefore, it follows: 

Lemma 3.7. For any q G M , let A(q) be irreducible. Then 1 is the maximum eigen- 
value and all the other eigenvalues A(q) G C are strictly smaller in absolute value of 
1, i.e. |A(q)| < 1. The eigenspace associate to the eigenvalue 1 is one dimensional. 
Furthermore, there exists an eigenvector rj G M N associated to 1 with strictly positive 
entries, i.e. A(q)r] = r\ and r\i > for i — 1, ..,N. Moreover, rji = -7= for i = 1, .., N 
and there are no other positive eigenvectors except multiples off). 

Proof Because for any q G M Nd , \\A (q) ||oo < rnax^i,..^ YujLi a «,i (q) = 1> we nave that 
the maximum eigenvalue is 1 and any other eigenvalue A(q) G C is strictly smaller in 
absolute value of 1. By Perron Frobenius Theorem the maximum eigenvalue is simple 
and the associated positive eigenvector is 77 with 77, = -7= for % = 1, .., N. □ 

Lemma 3.8. For any q G Mr, so that A(q) zs irreducible, let A(q) be the matrix 
as in A3. 15]) . We have that 1 G E(A(q)) zs i/ie maximum eigenvalue. The associated 
eigenspace is the d-dimensional manifold {p G R : p% = v, i = l,..,N;v G M. d }. All 
the other eigenvalues A(q) G E(A(q)) are such that |A(q)| < 1. 

Proof. It is an immediate consequence of (I3.24p and Lemma 13.71 □ 

We have finally the following result. 

Theorem 3.9. For any q G Mr , so that A(q) is irreducible, let C(q) &e defined in 
A3. 13]) . We have that G E(C(q)). T/ie (AT + l)d dimensional manifold X defined in 
A3. S\) is the eigenspace associated to the eigenvalue 0. All the other eigenvalues o/C(q) 
have real part strictly negative. 

Proof. From (13 . 1 9[) and Lemma [3TB1 we deduce that G E(C(q)) and all other eigenvalues 
have real part strictly negative. It is immediate to see that the algebraic multiplicity 
of is Nd + d. The (N + 1) d-dimensional manifold X defined in f)3.2p is the associated 
eigenspace. Namely, if w G X then C (q) w G X. From this one deduces that X is an 
eigenspace for the matrix C(q). Moreover, since the kernel of C 2 (q) is X, we get that X 
is the eigenspace associated to the eigenvalue 0. □ 

We denote by a(q) the spectral gap of the matrix C(q), that is 

a(q) := min {|Re(A(q))| : A(q) G E(C(q)), Re(A(q)) < 0} . (3.25) 

Let q G M Nd so that A(q) is irreducible. By Theorem 13.91 X is the eigenspace associated 
to the eigenvalue of C(q) for any q. We can therefore decompose M 2Nd as follows: 

M 2Nd = W(q)®l (3.26) 
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in such a way that W(q) and I are eigenspaces of C (q) and denote by n(q) the projection 
operator 

II(q) : R 2Nd -> W{q) . (3.27) 

3.3. Asymptotic Analysis. Let w° = (q°,p°) G Z be such that A(q°) is irreducible 
and let us set, for any r > and e > 0, 

B(r, e, w°) := {w = (q, p) G R 2Nd : |q - q°| < r ; |p - p°| < e} . (3.28) 

Denote by r the biggest value of r such that, for any w = (q, p) G B(r , e, w°), A(q) is 
still irreducible and «(q) > |«(q )- We set 

B(r„, e, w°) := ju; = (q, p) G S(r„, e, w°) : «(q) > ^a(q°) } . (3.29) 

The existence of r is granted since, by assumption, A(q°) is irreducible and Ur is 
smooth. To apply the spectral results obtained for C(q), (q fixed), to the nonlinear 
system f)3.12p we write 

C(q(f))=<7(q(0))+r(q(f)) , (3.30) 

where 
and 

fl(q(f)) := A(q(t)) - A(q(0)) . (3.32) 
Next we estimate the norm of U(q(£)). 

Lemma 3.10. Let (q(£),p(£)) fre i/ie solution of A3.12\) starting from the initial data 
(q°,p°). Then, setting 

Mq(*)) :=Oij(q(*)) -Oij(q°) * 5 J = , ( 3 - 33 ) 

||S(q(t))||<2iV^f^Mi sup |-(g? - $ + *(*)- . (3.34) 

We defer the proof of this result to the appendix. 

Theorem 3.11. Let w° = (q°,p°) G X and assume that A(q°) irreducible. There 
exist three positive constants r = r (u>o), T = T (w°) and €q = eo (iV, w°) such that, for 
any initial datum w 1 G B(r ,e ,w°) C R 2Nd , 

dist (wfow 1 )^) < e e-*^ , (3.35) 
where w(t,w r ) is the solution at time t of AS.lfy) starting from w 1 . 
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Proof. For any s > 0, we define 

Q(s, w°) := {w = (q, p) G R 2Nd : | [1^ - 0\ (q - q°) | < s} , (3.36) 

where f2 is the operator defined in (13. 4p . Denote by so the largest value of s such that, 
for any w = (q, p) G Q(s , iu°), ^4(q) is still irreducible and a(q) > ^a(q°). Such a value 
s exists since A(q°) is irreducible and Ur is smooth. Let us set 

Q(s ,w°) := = (q,p) G Q(s ,w°) : a(q) > ^«(q°)} • (3.37) 

We have that 

B{r , e, w°) C Q{s , w°) , Ve > . (3.38) 
Namely we have that sq > tq since requirement (13.291) is stronger than fl3.37j) and 

Ifl^-OKq-q )! < |q-q°| <r . (3.39) 

Let w(t, w 1 ) = (q(t, w 1 ), p(t, w 1 )) be the solution of system (I3.12p starting from an 
initial datum w l G B(r Q ,e,w°) and let t*(wi) > be the first exit time of w(t, w 1 ) from 
Q(s , w°). If w(t, w 1 ) G Q(s , w°) for all t > 0, then we set t*(w r ) = oo. Next we analyze 
the solution for t < t*(w l ) and we will show that t*(w 1 ) = oo for any initial datum 
w 1 G B(r ,e,w°), provided that e in ( I3.28P is suitably chosen. Let us define 

^-n^^M^ 1 ) , (3.40) 

X(t) ■= (h Nd - n(q(t, w 1 )) w(t, w 1 ) , t< t^w 1 ) . (3.41) 

By construction x(t) G X, £(t) G iy(q(t, w; 1 )). We then have 

= (^W))) w (*> + n(q(t))^(t, w 1 ) (3.42) 



d H(q(t)) ) «,(*, «,!) + II(q(t))C(q(t))u;(f , W 1 ) 



dt 

Taking into account that w(t, w 1 ) — + x(t) we get 

Jt m = (it n(q(t)) ) m + (l n(q(t)) ) + nfaWM^))^) • ( 3 - 43 ) 



d 

111U11 ll^ll^fcJJ, we nave 

C(q(£)) and II(q(t)) commute, we obtain 



Since for any given w G X, by the definition Il(q(t)), we have ^n(q(t))w = and 



^) = (| n («l(*))) «*) + C(q(*))eW • (3.44) 

Setting 

C(q(t)) = C(q(0)) + r(q(t)) , (3.45) 
where T(q(t)) is defined in (13.31 j) . we get 
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By the formula of variation of constants: 

m = e c(q(o))t£ (0) + j\cwow-s) j^ n (q(s))) e(a) + r(q(a))^(a)}d«. (3.47) 

Performing the exponential of the matrix C(q(0)) one needs to take into account that, 
because of the possible presence of Jordan blocks, powers of t might appear. We control 
such terms paying e - ^^ ))* and multiplying the remaining exponential by a constant 
D(C(q(0))) which depends only on C(q(0)). Since q(0) G B ro (q°) , which is a compact 
set in M. Nd , we denote by D := sup qgBrf) ( q o) Z?(C(q)), which depends only on C(q°) and 
r . Therefore, we get 

< C o e-^ qW>t |e(0)|+^ jTe-i^t*--') 1 1 Q^II(q( a ))) £(s)| + |r(q(a))£( a )|}<fo 

(3.48) 

Next we estimate ||^n(q(t))|| . Let II(q(t)) = {^-(q^))^ . =1 ^ , we then have 

^n(q(i)) = (V q(t) 7r f>i (q(t)) • p(/)) ; ; , v (3.49) 

= (V q ( t )7r^-(q(t)) ■ [p(t) - fip(t)]). . =lj _ i>JV . 

The last equality holds since, by flBTTTj) . II(q(t)) = n(q(0)),Vi G M, when q(i) is the 
evolution given by the flow on the invariant manifold, i.e. when p(£) = Qp(t). We get 
by Corollary 13.51 



< sup sup |V g 7Tjj( q )| |p(t) — Op(t)| (3.50) 

{qeK^d : «j=(q,p)eQ(s ,™ )} i,i=l,..,JV 



<£>'(s Q )e, VtGfO,^ 1 )) , 
where D'(so) > 0. Furthermore, by f 1 3 . 3 1 j) and Lemma [3.101 we have 

||r(q(t))|| = \\B(q(t))\\ < 2-^— sup \VU(x)\ max \-{q} - q\) + - q k (t)\ 

UR (0) a-gjjd l<i,k<N 

(3.51) 

and, by Theorem 13.4} for i, k — 1, .., N, 

\qi(t)-q k (t)-(ql-ql)\= [ \ Pi (s') - p k (s')\ ds' (3.52) 

Jo 

= [ \ Pi (s')-p i +pl-p k (s')\ds'<2et, 
Jo 

where p° = p° k since (q°,p°) G X. Thus, setting D 1 := 2 ^*gg^ , Vt G [O,^ 1 )) we 
obtain 

|£(t)| < D e-3 a ^>*|e(0)|+Doe /* e -Hi(°))('- s ) {[£>' ( s „) + 2jVL> lS ] . (3.53) 

Jo 
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Given K > max{_D' (s ) , 2Di}, take T G (0, ^(w 1 )). A suitable choice of T will be done 
later. Then, 

|£(t)| < Doe-5 a ( q (°W t |^(0)| + eDoA'{l + TJV} /V^ (q(0))( '- s) |£(s)|ds, Vt G [0,T]. 

(3.54) 

By the Gronwall's inequality, which we recall in appendix (see Lemma [5. ip . we get 

\£(t)\ < L>o|e(0)|e-*^ a(q(o)) - e5 ] < D |£(0)| e -*[s a(q °M , Vt G [0,T] , (3.55) 
where we made use of ( 13. 3 7ft and set 5 := D K {1 + AT} . Let us choose e such that 

^-a(q°)>eD K{l + NT}. (3.56) 
lb 

Then, 

\m\ < D \t(0)\e~^ a{cl0) , Vt G [0,T] . (3.57) 

Since 

dist (w(t,w l ),l) = inf w 1 ) - w\ = inf |x(t) + £(£) - w| (3.58) 
wex wax 

<|£(t)| + ini Mx(t)-w| = |e(*)l, 

we have 

dist (wfow 1 ),^) < A)|£(0)|e-^ a(qO) , Vt G [0,T] . (3.59) 
Then, recalling that dist (w(t,w l ),T) = |p(t) — Qp(t)\, and, since for w 1 G B(r ,e,w°), 

|£(0)| < D(s ) dist (V,X) < .0(50)6 , (3.60) 

we have 

|p(t) - fip(t)| < OoO(s )ee-^ a(q0)t < e(D Q D(s ) V l) e -^ Q(q0) * , Vt G [0,T] . (3.61) 
From this we get 

\[I Nd - Q] (q(t) - q°) I < \[I Nd - 0\ (q(0) - q°) | + / | [1^ - fi] p(s)\ds (3.62) 

< | q(0 ) _ q 0| + (L> L>( So ) V l)^^ 1 " e^^) 

a(q°) 

1 6 

<|q(0)-q | + (A>O( S0 )Vl))^ e . 

Let us choose e such that 

16 1 

r + (D D(s ) V 1 W^e < -s Q , (3.63) 
and denote this chosen value by e\. Now we first choose T such that 

(D D(s ) V l)e-^°) T = i A ^1- (3.64) 
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and denote this chosen value by T , then we choose e 2 in such a way that ( I3.56P holds 
with T replaced by T . We then set 

e := mm{e 1 ,e 2 } ■ (3.65) 

Notice that, by f l3T63|) . 
so, since a (q°) < 1, 

eo < h < . (3.67) 

We remark that the choice of To depends only on w° G X, while the choice of eo depends 
on w° G X, on T and N. Therefore, at time T we have 

\[I Nd -n](q(T )-q°)\<Ls (3.68) 



dist (w(T , w x ),X) = |p(T ) - fip(T )| < | . (3.69) 



and 



We can then repeat the previous argument for the solution of the system (12.1 ft starting at 
time T from the initial datum (q(T ), p(T )). We need to recall that a(q(T )) > |a(q°). 
In a similar way we can show that that for t G [To, 2To], 

|p(f)-Op(f)| < D \aT )\e-^-^ a ^ , W G [T ,2T ] . (3.70) 
Therefore, by ( I3.64p . we have 

dist M2T XTo)),X) = |p(2T ) - fip(2T )| < e ° < § , (3.71) 

2 2 (A)T'( s o) V 1) 2 2 

and, by f l3T70|) . 

|pW d - 0\ (q(t) - q°) | < Ifljvd - (q(T ) - q°) | + / | [Ijvd - fi]p(s)| tZs (3.72) 

■/To 

< - fl] (q(T ) - q°) | + D \aT )\ f e^-^^ds 

1 1 

~ 2 S ° + 4 S ° ' 

the last inequality being a consequence of f !3.63|) and (I3.64p . Thus, at time T x = 2T 

|[E OT -n](q(r 1 )-q )|<i So + ^. (3.73) 
Hence, we have that (q(Ti), p(T)) G Q(s , w°). Iterating this procedure m times we get 
dist (w(T m ,w l ),Z) = |p(T m ) - fip(T m )| < ^ , (3.74) 
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and 

| [I Nd - Q] (q(T m ) - q°) | < s ]T • (3.75) 

fc=0 

Since Efc>i 2^" = I we obtain the thesis of the theorem. □ 

4. Kinetic limit: Boltzmann-Vlasov equation 

We would like to study system (12. 1 p when the number of particles N goes to infinity 
and to derive the kinetic equation for the density ft{x, v) of particles at x with velocity 
v at time t. Let (qj(t),pj(t)), j = 1,..,N be the solution of the system (12. ip for some 
initial datum (g°,p°), \\p°A\ < l,j — 1, -,N and 



v 

.N 

I ll.l . II I I - : 

N 



(4.3) 



//^ (dx, dv) = — S(qj(t) — x)5(pj(t) — v)dxdv , (4.1) 

the empirical measure. By Lemma 13.11 and Remark 13.21 fi^(dx, dv) has support on 
M. d x Bi. Given a smooth function g on R d x B±, we denote by 

H?(g) = [ g(x,v)fj?{dx,dv) . (4.2) 

JR d xB 1 

Writing the second equation of (12. ip in term of we get 

dpijt) _ Ef=i uMf) - qM (pM - Pi (g) 
* -<&(*)) 

_ J RdxBl - j/) (m - Pi (*)) (M?(dy, du) 

J RdxBl U R (qi{t) - y)fi?{dy, du) 
=: M( qi (t), Pi (t),Li?) . 
Therefore, the evolution of iif is given by 

= ^ (v . Vxg) + tf {M ( . ; . ^ j . (44) 

In the equation (14.41) . iV is fixed. To study the limit as iV — > oo we assume that at t = 

1 r 

/i^ (<ix, rff) = — 5(g° — x)5(p® — v)dxdv ==> fi (dx, dv) , (4.5) 

where the convergence is weakly as measures, i.e., for every bounded and continuous 
function g, 

lim / g(x,v)no (dx,dv) = \ g(x,v)fj, (dx, dv) . 

N^ooJ J 

We want to show that if ( 14. 5 p holds at time t = 0, then 

fxf(dx,dv) Ht(dx,dv) , (4.6) 
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where (x t is the measure solution of the following equation 

= ih(v ■ V x g) + fi t (M (-, -, ■ V v g) , (4.7) 

which is the formal limit of (14.41) . Furthermore, if 

/i (dx,dv) — f (x,v)dxdv , (4.8) 

i.e. the initial measure is absolutely continuous with respect to Lebesgue measure, then 
we want to show that 

Ht(dx,dv) = f t (x,v)dxdv , (4.9) 
where f t is the weak solution of 

^f t (x, v)+v V x f t (x, v)+V v - [M(x, v, f t )f t (x, v)] = 0. (4.10) 

If one assumes that fit(dx, dv) = ft(x, v)dxdv, equation (I4.10p is obtained integrating by 
parts in (14.71) 

d f 

—fi t (dx,dv)g(x,v)= v •V x g(x,v)fi t (dx,dv) (4.11) 

xB, JR d xB 1 



+ / M(x,v,/i t ) ■ [V v g(x,v)]/i t (dx,dv) 

'xBi 



Therefore, 



d f 

—f t (x,v)g(x,v)dxdv = - v ■ V x ft(x, v)g(x, vjdxdv (4.12) 

g(x, v)V v ■ [(M(x, v, f t )) f t (x, v)} dxdv , 

where, for F(v) E R d , we denote by V v ■ F(v) = E?= i = div(F(v)). This 

holds for any test function g(-, •); therefore equation (14.101) holds. This is a sort of 
Boltzmann-Vlasov equation where the collision kernel is replaced by M(x,v, f t ). When 
M(x, v, ft) = the equation ( 14.10p becomes linear and describes the motion of particles 
moving independently from each other with the same velocity. 

We prove rigorously the previous argument under some assumptions over the inter- 
action Ur and the configurations space. We first show in Theorem 14.51 the existence in 
a suitable space of measure of the evolution (14. 7p . Then we show in Corollary 14.71 that 
if ( 14. 5 p holds, then ( 14. 6 j) holds and it is the measure solution of ( 14. 7p . Furthermore, in 
Theorem 14.81 we show that if ( 14.81) holds then (14.91) holds and ft is the weak solution 
of (14. lUp . Requiring more regularity to the initial distribution f and to the interaction 
Ur we proved that f t is the strong solution of (I4.10p . 

The results are shown adapting to our context the method reported in Spohn's book 
(1991) [S] (see also Neunzert (1984) [N] and Dobrushin (1979) [D]) and some classical 
tools of dynamical systems. 
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In Subsection 14 . 1 1 we present some qualitative properties of the solution of (14. TJ) . More 
accurate analysis of the qualitative behavior of this solution, in special way its long time 
behavior, is behind the aim of this paper. 

We start defining the space of measures and the metric we will be using. We will 
consider measures with finite total mass on (X x B\, B (X x B\)) where the symbol X 
stands either for M. d or for the torus of linear size D > 0, Td and B (X x Bi) is the Borel 
aalgebra on X x B\. We will denote by M. the space of such measures and, without loss 
of generality, we will take into account only probability measures. Notice that we will be 
using the same notations either to denote the space of probability measure on Td x B\ or 
the space of probability measure on R d x 5], unless we will have the need to distinguish 
between the two configuration space and therefore use the notation M.{X x Bi). The 
bounded Lipschitz distance^, d bC , between two measures fi and v in M. is given by 



<4c(/i, v) = sup 



g(x,v)/i(dx, dv) — / g(x,v)p(dx,dv) 



(4.13) 



where 



V:=->g\g:XxB 1 ^[0 } l] } \g(x, v) - g(y, w)\ < J\v - w\ z + \x - y\ z \ . (4.14) 



For any smooth function g we consider the evolution (14. 7ft in M.. To prove existence 
and uniqueness of the solution of equation (14. 7J , an important assumption we make is 
that M(-,-,z/) is Lipschitz continuous in x and v for any v G M. We are then forced 
to consider only those interactions Ur for which not only M(-, -, v) can be well defined 
for all (x,v) G X x B\ and for any v 6 M, but is also Lipschitz continuous in x and 
v for any v G M.. As a matter of fact, it is not clear even how to define the quantity 
M(x,f,/if r ) for (x,v) G X x B\ regardless of how smooth one can take the interaction 
Ur. The problem is caused by the presence of the denominator in (14.31) . When x is in 
the support of fif the denominator is always positive, actually bigger than and then 
M(x,v, fif) is well defined. When x is not in the support of it might happen that 
J R d xBl Ur(x — y) (dy , dv) = 0, therefore it is not clear how to define M(x,v, To 
overcome this problem we consider two classes of interaction Ur. The first one is the 
collection of smooth interactions Ur such that (U r* v)(x) > for all x G X. In this case 
we define 

,w v I J XxBl UR(x-y)uv(dy,du) \ 

M(x,v,u)=\— — — \ - v-. A(x, v) —v . 4.15 

V Jxxb, Ur{x -y)u{dy,du) J 



2 The bounded Lipschitz distance is identical to the Kantorich-Rubinstein (Vaserstein) distance. The 
metric dbc generates the weak *-topology on M. For a sequence fi N , 

lim d bc {fJ, N , v) = 

N— ¥oo 

is equivalent to 



\im j g(x,v)ii (dx,dv) = j g(x,v)v(dx,dv) 



for all bounded and continuous functions g. 
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The second one is the class of smooth interactions Ur with compact support. In this 
case, we fix e > and define 

M (x,v,v) := . . (4.16) 

JxxB 1 u R( x -y) iy ( d y^ du ) + e 

Hence, in the case Ur has compact support, we modify the interaction term in such a 
way that when f XxB Ur(x — y)u(dy, du) = then M e (x, v, v) = 0, when f XxB Ur(x — 
y)u(dy, du) > e then M e (x, v, v) = M(x, v, u)+0(e), when e > J XxB UR(x—y)b , (dy, du) > 
then M e (x, v, v) is a large perturbation of M(x, v, v). Notice that if one considers the 
system (12. II) replacing Ur with the modified interaction, the one in ( 14.161) . then Lemma 
13. II continues to hold. This is the reason why we kept the velocity in the ball B\. 
It is easy to see that for any measure v on IR d x B\ we have 

sup \M(x,v,u)\ < 2 , sup \M e (x,v,v)\ < 2 . (4.17) 

(x,v)eR d xB! (x,v)£R d xB! 

The Lipschitz continuity of M(-,-,v) with respect to v follows from the linearity of 
M(-, •, v) as a function of v. The Lipschitz continuity of M(-, •, v) with respect to x does 
not hold in general even if one takes a smooth interaction Ur. The problem is created 
by the presence of the denominator in M(-,v, v). The denominator of the derivative in 
Xi,i = l,..,d, of any component of the vector A(-,u) might be very small while the 
numerator, because of the presence of the derivative of Ur, might be not of the same 
order. Hence, one can certainly control the gradient of A(-, v) if the gradient of Ur is of 
the same order of Ur. We have then the following lemma. 

Lemma 4.1. Let v G Ai and let K be a positive constant such that sup xgK d |Vlog[/#(a;)| < 
K. Then 

\Ai(x,u) - Ai{z,v)\ < L\x-z\ , x £ X,z 6 X, i = 1, ..,d, L = 2K . (4.18) 
Proof. Vi = 1, .., d, we have 

f XxBl UiV x U R (x - y)v(dy,du)(U R *v)(x) - f XxBi U R (x - y)uiis(dy,du)V x (UR*v)(x) 
V x Aj(x, v) — — ■ 

[(f/ R ^z/)(x)] 2 

(4.19) 

/(xxfii) 2 K^f > du)u(dy', du')ui [V x U R (x - y)U R {x - y') - C/r(x - y)V x U R (x - y')] 
= 1(Ur*v)(xW ' 
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Taking into account that \u\ < 1 we have 

LxxBiV v(dy,du)v(dy',du') \V x U R (x - y)U R (x - y') - U R {x - y)V x U R (x - y') 



\V x Ai{x,v)\ < 



< 2 



\(U R *v){x)Y 

■f(xxBi) 2 u ( d V^ du)v(dy', du') \V x U R (x - y)\ U R (x - y') 



[(U R *v)(xW 

/(XxBi)2 v(dy,du)v(dy',du') \V x U R (x-y)\ ^=^U R (x - y') 



[{U R *v){x)f 



(4.20) 



<2sup l V ^-f)l< 2 ^ 
y U R {x-y) 



□ 

Remark 4.2. The assumption about the boundedness of sup xeR d |V logU R (x)\ is quite 
strong. As a matter of fact, an interaction U R verifying this assumption should decay 

]x\ 

for \x\ large as e~~R . Interactions with compact support do not satisfy this assumption 



as well as the interaction M d 3 x i — > U R (x) 



T e 2R 2 e 



considered next. 



Lemma 4.3. Let To be a d dimensional torus of linear size D > 0, 

1 



3 x i — > Ur(x) 



-e 2R 1 g 



(4.21) 



(2ttR 2 )% 

and V R be the corresponding interaction defined on Tb through the periodization of U R 

M. d 3 x i — > V R {x) = U r( x + nD ) e . (4.22) 



For v G M. and A(-, v) as defined in h4-15\ ), with U R replaced by V R we have 



\Ai(x, v) — Ai(z, v)\ < L\x — z\ 



x,z<eTd, i = l,..,d, L = . (4.23) 



Proof. Let us write 

\Ai(x,fi) - Ai(z,(jL)\ 



f 1 d 

/ ds—AAsx + (1 — s)z, /i) 
Jo ds 



< sup 

se[o,i] 



— Ai(sx + (1 - s)z,/j.) 
ds 



i = l,..,d, 



(4.24) 
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and set x = sx + (1 — s)z and (V R * fi)(x) = J rDXBl Vr(% — y)f^(dy, du), we obtain 
d., , (x - z) !t d xb 1 ~ {X r~ V) v r( x o - y) u iKdy, du)(V R * /i)(x ) 

Ts Mx °' " ] = ~R PWOM? (425) 

(x - z) fr DX B 1 v r( x o - v'KKdy', du') J TdxBi ~ {x R ~ y) V R (x - y)^{dy, du) 



R l(V R *i2)(x )]i 
(x — z) 



R 

where we denote by 



(4.26) 

Ci(x ,R,fx) , (4.27) 



d{x Q ,R,ii) := \ x (4.28) 

[Xn — y ) 

-V R (x - y)uifi(dy,du)(V R -k fi)(x ) 



Vr(zo - y^u'^y', du) ! X ° V V R (x ) [i(dy , du) 

T D xBi JTdxB! K 

Finally, rewriting 

C (x R ) = / {j ^ 1 Vb{xq - y') Vr(x - y) J u'^jdy', du')n(dy, du) ^ 
JV R (x -y')V R (x -y) J fi t (dy',du')jj,(dy,du) 

and recalling that \u\ < 1, we obtain 

\C{x R fi)\ < ^ TdxBi ^ Vr ^ x ° Z y '} VR ( x ° ~ y } U d xb 1 du')fi(dy, du) ^ d 

°' ' ~ /r D xBi ^(^o - y')VR(a;o - 2/) J TDXBl fJ,{dy',du')fLt{dy,du) ~ R 

(4.30) 

since in the torus \y' — y\ < D and the result follows by (I4.25p . □ 

M € (-,v, v) is easily seen to be Lipschitz continuous in X. In fact we have: 

Lemma 4.4. Let v e Ai,e > 0,U R (-) a smooth interaction whose support contained in 
a ball of radius R such that sxvp xeB \VU R (x)\ < 1 and M £ (-,t>, v) as in Then 
M e (-,v, v) is Lipschitz continuous in X. 

\M*(x,v,v)-M*(y,v,v)\<L\x-y\, x,y e X,i = 1, ., d, L = ~. (4.31) 

To prove the existence of the solution of (14. 7p we prescribe a curve t — » fi t G M. 
weakly continuous in t and we consider the following non-autonomous system of ordinary 
differential equations: 

i<t) = v(t) 

iv(t) = M(x(t),v(t)^ t ) ■ [ ^ 6Z) 



DYNAMICS AND KINETIC LIMIT FOR VICSEK-TYPE PARTICLES 



21 



Under the assumption that M(-, -,//t) is Lipschitz continuous in X x B\ there exists an 
unique global solution of (I4.32p for any given initial datum. The corresponding time 
dependent two parameters flow is denoted by Tt jS [//.]. Under this time dependent flow 
any initial measure evolves as 

v t = v °T ,M , (4.33) 
where To 5 t[/i.] is the push forward of the measure z/o. For any test function g we have 
that 

"t(j) = ^o3t,oH) , (4-34) 
where T tj0 [/i.] is the pull back acting over the test functions. By the existence and 
uniqueness of the solution of ( I4.32p for any initial datum, the inverse flow (T^J/x.]) -1 is 
well defined. The equation for the evolution of z/ t , easily derived, is 

d(M9)) 9(u (goT t0 \p.})) 

— Q- f — = ^ = M( v ^x9) ° T ti o[n.\) + u ((M(x, v, ii t )V v g) o T tfl [^.\) 

(4.35) 

= v t (v ■ ^x9) + "t(M(x, v, fx t ) ■ V v g) . 

One immediately realizes that proving the existence and uniqueness of the solution of 
(14. 7p is equivalent to prove the existence of a fixed point for the time dependent flow 
/ij = |io o T 0i t [/-£■]• This is the content of the next theorem. 

Theorem 4.5. Let Ur be as in Lemma \4-l\ or as in Lemma \4-^\ and let M(-,-,z/) be 
defined as in Hi4-15\ ) for any v G -M(7b xB^. The equation (TJT^j has an unique solution 
in the space A4(Td xBi) if fi G M.{Td xBi). Furthermore, take two solutions of M-ty , 
fit starting at fiQ = fi and v t starting at uq = v then in the bounded Lipschitz distance 

d b c(v t , Vt) < e ct d bC (fi, v) , (4.36) 

where c is a constant which depends on the Lipschitz constant of M(-,-,u) and on 
inf xe7 - D U R {x) =: a > OH 

The proof is obtained adapting the method explained in Chapter 5 of [S] to our context. 
The main difference between the case considered here and the one presented in Section 5 
of [S] is that, in our case, the dependence of M(-, •, v) from v is not linear. To overcome 
this problem we assume that the interaction Ur is such that inf xg 7- D Ur(x) = a > 0. This 
is the case for interactions as in Lemma 14.11 but defined on Tb x B\ via periodization 
or as in Lemma 14.31 To facilitate the reader we report the proof of Theorem 14.51 in the 
Appendix. 



Remark 4.6. Theorem |^.<5| does not hold in M. d x B\ when Ur satisfies Lemma 4-1 
Although in this case Ur is globally Lipschitz continuous in M. d , we are not able to show 
that M(x,v, •) when x G M. d , is Lipschitz continuous with respect to v G M. in the d^c 
metric. The theorem applies with obvious modification if we take an interaction Ur with 
compact support and define M € as in ( [^.i6] ). In such case the theorem holds either for 



3 Notice that for Ur as in Lemma [4.11 or as in Lemma [4.31 (Ur * v)(x) > a for all x e Tb, where 
a = w£ xe r D Ur[x) > 0. 
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the system defined on Tb x Bi or on ¥L d xBi. The constant c in the statement of Theorem 
\4-5\ will then depend on e, the lower bound of the denominator of M e . 

Corollary 4.7. Let ji G M. and (cf N ,p N ) G (To x Bi) N be a sequences of particle 
configurations, so that 

lim d bC (fJ. N , fj) = , (4.37) 

AT— s>oo 

where fi N is the empirical measure. Let (qjv(t),pjv(^)) £ (7d x B{) be the solution of 
112. 1\) with initial datum {c[%,p%) and with Ur chosen as in Lemma \4-1\ in Td x B\ or 
as in Lemma \4^3\ Then there exists |ii6M such that 

lim 4c(/if,/i t ) = , (4.38) 

N— >oo 

and \i t solves equation 

The proof is immediate from (I4.36P . The validity of Corollary 14.71 for smooth compact 
potential Ur and for the local mean velocity increment M e is immediate as well. 

Theorem 4.8. Let M(-, •, fx) be as in d4-15\} and assume that M(-, •, jj) G C 1 (X x Bi) for 
li G M.. If fio(dx,dv) = fo(x,v)dxdv, then ii t (dx,dv) = f t (x,v)dxdv and f t is the weak 
solution of MJVj) . Furthermore, if f G C k {X x B x ), k > 1, and M{-, •,//) G C k {X x B t ) 
for LieM, then f t G C k (X x B ± ). 

Proof. We start showing that for any given weakly continuous curve t — > fi t G M, 
if vo(dx,dv) = q (x,v)dxdv, i.e. absolutely continuous with respect to the Lebesgue 
measure, then u t (dx,dv) = q t (x,v)dxdv, where 

d 

— q t (x, v) + V x q t (x, v)-v + V v - [M(x, v, Li t )q t (x, v)\ = , (4.39) 

and, if q G C k (X x Bi) and M(-,v,fi) G C k (X x B ± ) for any /i G M, then g t G 
C fe (X x B\). Note that (14.391) corresponds to a linearization of (14.101) since M(x,v,fj, t ) 
does not depend on q. once /if is given. In Theorem 14.51 we proved that the fixed point 
equation fi t = yUo ° ToAt 1 -] holds. Therefore, by this result and the validity of (14.391) . one 
immediately obtains that Lit has density and the thesis of the theorem is proven. We are 
then left with the proof of (14.391) . Let us set w — (x, v) G X x B\. For any test function 
g we obtain 

Mg) = T 0it \fi.](g) = vo(g° T t>0 \fi.]) (4.40) 



u (dw)(g o T tfi [^])(w) = / q (w) (g o T tfi [fU-}) (w)dw 
q (w) o (T t)0 [ii.])~ 1 J(w,iit)g(w)dw 

'JxBi 

where J{w, fi t ) = Det [9.(T ti0 )[//.]) _1 (w)] is the Jacobian of the flow (T t0 [/i.]) _1 computed 
in w. Since the divergence of the vector field (v(s), M(x, v, /i a )) is given by 



E 



dvj M(x,v,fx s ) 
dxj dvj 



-d , (4.41) 
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by Liouville Theorem (see |Alj or |A2j ) for any weakly continuous curve t -y /J t G M, 
we have 

Det [d.(T tt0 )\pt.])] (w) = e~ dt , E X x B, , (4.42) 

hence, 

J(w, fi t ) = e dt , VtoGlx^. (4.43) 
Then, from (jOD|) and (jQ5|) . we obtain 

z/ t (#) = / e d *(g o(Ti, [/i.])- 1 )H(?H^= / gtHsHdw, (4.44) 

JXxB 1 JXxB 1 

where we denote by 

q t {w):^e dt (q o{T t ^])- 1 ){w) . (4.45) 

Notice that qt{w) is weakly continuous in time, since (i. is weakly continuous. Fur- 
thermore, if, for k > l,M(-,-,/i) G C k (X x Bi),fi E M and g e C k (X x 5^, then 
gt G C fe (X x Si). Writing e~ dt (q t o T t;0 \p.])(w) = qo(w) and differentiating with respect 
to t we get 

^ (e~% o (T ti0 [M.])H) = -^- d '(% o T t , [ M .])H+ (4-46) 

+e" d *^(g t o r ti0 [/i.])(w) + e-*V s ( ft o r t)0 H)H • u ° T tj0 H+ 
+e-*V w ((ft o T t)0 [//.])H) ■ (M(; ; /it) o T t) oN) W = . 

Multiplying both members of the previous identity for e dt and applying to (T| ) o[/i.])~ 1 iu 
we obtain 

d 

- dq t (x, v) + -7^qt(x, v) + S/ X q t (x, v ) • v + V v q t (x, v) ■ M(x, v, fi t ) = . (4.47) 

Notice that this last equation is linear in q. since fi. is given. Therefore, the equation for 
ft is 

^f t (x, v) - df t (x, v) + VJ t (x, v)-v + M(x, v, f t ) ■ V v f t (x, v) = (4.48) 

which corresponds to ( 14. 10ft . □ 

Remark 4.9. Theorem \4-S\ can be proven when Ur is a smooth function with compact 
support and the local mean velocity increment is M e for e > as defined in ^.lo] ). In 
this case the theorem holds either in W d x B\ or in Tb X B\. The only difference in the 
computations done in Theorem \4-8\ is that 

V-M e (x,^ s ) = -d ,^**^*] = -dh s (x) = -dh{n s ){x) . (4.49) 

Then, for any weakly continuous curve t — > fit G M. and any w G X x B\, by Liouville 
Theorem, we have 

Det [<9.(T tj0 [/i.])] H = e -d.^ds{h s oT s ^.])( w ) ^ (4 _ 5Q) 
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therefore 

J(w,nt) = Det [atCT^oH)- 1 ]] = e d ^ d <hsoT s , t[fl .])( w ) _ ^ ^ 

Then, from h4-40[ ) and h4-51\ ), we obtain 

Pt{g) = [ e dftMhs°TMM ( go o (T^.])- 1 ) (w)g(w)dw (4.52) 

JTdxBi 

q t {w)g(w)dw , 

T D xBi 

where we denoted by 

q t (w) := e if!Mhs°TsMM ^ (T^.])" 1 ) (tu) . (4.53) 

Notice that qt(w) is weakly continuous in time, since /x. zs weakly continuous. Further- 
more, ifM(-,-,fi) E C k {X x B 1 ),fx E M andq E C k {X x B 1 ), then q t E C k {X x B x ). 
Writing 

e - d IoMhsoT M ( W ) {qt T t M)(w) = q (w) (4.54) 
and differentiating with respect to t we get 

+ e -d/>(^o(T s , [M.])HfL(^^^^ 

+ e-^o ^(«») V , o T ti0 [//.])H) ■ (M(-, ^) o T, [/i.])( W ) = . 

Multiplying by e d fo ds ( hs °( Ts -°^( w ) and applying to {T t $[[i])~ l w we obtain 

d 

- dh(^ t )(x)q t (x,v) + —q t (x,v) + V x q t {x,v) ■ v + V v q t {x,v) ■ M(x,v,fx t ) = 0, (4.55) 

which is linear in q. since /i. is given. Therefore the equation for f t is 

^f t (x, v) - dh(f t )(x)f t (x, v) + V x f t (x, v)-v + M(x, v, f t ) ■ V v f t (x, v) = , (4.56) 

which corresponds to l[4-10\ ) taking into account the definition of h in ^-4^ - 
4.1. Qualitative behaviour of the solution of (14. 7h . 

Lemma 4.10. Let t — > fi t £ M- be the solution of fl^. 7| ) with initial datum hq. We have 

Ht(x) = hq(x) + / n s (v)ds , (4.57) 
Jo 

\v | 2 fit{dx, dv ) < / \v\ 2 fio(dx,dv) . (4.58) 

XxBi JXxB 1 
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Proof. By (14. 7p we have 
d_ 
dt 



Xi/j, t (dx, dv) = [i t (v. 



XxB 1 



which imply (14.571) . To obtain (I4.58p . again from equation (14. 7p we get 



d 



\v\ 2 fi t (dx,dv) = ^MM 2 ) 



(4.59) 



(4.60) 



fh(v ■ V x \v\ ) + /i 4 (M(-, •, fx t ) ■ V v \v\ ) 
ij 



= 2^/i t (M i ( v ,// t H) <0 . 

i=l 

Namely, for i = 1, .., d, when Mj(-, ■, fit) ^ 0, we obtain 

fi t {dx, dv)Mi(x, v, nt)vi 

IxxBx U r( x ~ y) ( y i u i ~ v i) Mdy, du) 



fit(dx, dv) 



(4.61) 



XxBi 



XxB 1 



SxxBt u r( x - y)vt(dy, du) 

(A A \ (A A \ Ur ( X ~~ y^ Ui 

H t (dx,dv)n t (dy,du)- — . 

(XxBi)a JXxB^R^ -y)Mdy,du) JxxB 1 

by Schwartz inequality. 



Ht(dx, dv)v\ < , 

□ 
Any- 



Jensen inequality and f )4.58j) imply the boundedness of f XxBi Vfi t (dx,dv) 

way, we are not able to show more interesting estimates about the mean velocity 
JxxB l v Vt(dx,dv). 

Lemma 4.11. Take M as in fl^.i5| ). The equation is not time reversible, i.e. 

invariant under simultaneous reflection t — > — t and v — > —v. 

Proof. Let us set b t {x,v) := f_ t (x,—v). We have that j^b t (x,v) = —4rf- t (x,—v) and 
Sbtfx^ __ _ df- t (x,-v) ^ Therefore; by 



dvi 



— b t (x, v) = -— /_ t (x, ~ v ) = -v-V x f- t (x, -v)-df„ t (x, ~v)+^2 M i( x , -v, f-t)- 

i=l 



dt 
and 



dvi 
(4.62) 



M(x, -v,f_ t ) 



J K d xR d Ur(x -y)(u + v) f- t (y, u)dydu 

IxxB! U r( x ~ y)f-t(y,u)dydu 
Ixxb 1 U r( x ~ V) (~ u + v ) b t(yi ujdydu 



(4.63) 



Ixxb 1 u r( x ~ y)Hy^ u )dyd 



u 



-M{x,v,b t ) . 
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Thus, the equation for b t (x,v) is 

b t {x, v) + v V x b t {%, v) + db t {x, v)~22 Mi{x, v, b t ) — = 



dt 



which differs from equation (14.101) . 



8=1 



dvi 



□ 



Remark 4.12. Lemma 4-H also holds when Ur has compact support and M is replaced 
by M e . 



Let ft be the solution at time t of the equation (14.101) . We denote by H(f t ) the 
Boltzmann-Vlasov entropy 



ft(x, v) ln(f t (x, v))dxdv 



(4.64) 



XxB 1 



In the next lemma we show that H(f t ) is a decreasing function of the time. 
Lemma 4.13. Let f. be the solution of I{4-1U\ ) with M chosen as in ( [^.i<5| ), then 

-d . (4.65) 



Let f e be the solution of fl^.lOp with M replaced by M e chosen as in fl^.ifip , then 

d 



dt 



H(fl) = -d / ht{x)ft{x,v)dxdv 



(4.66) 



XxBi 



where, as in \4-49j )> ht 



Proof. We start showing (14.651) . The proof of (14.661) is similar and we will only outline 
the differences. 

±H(f t ) = - [ ^(x,v) [Hf t (x,v)) + l]dxdv (4.67) 
dt Jxxb 1 dt 

= / [In f t {x,v) + l][vVJ t {x,v) + V V ■ M{x,v,t)f t {x,v)]dxdv . 

JXxB! 



Integrating by part the last term in (14.671) we get 
d 



dt 



H(ft. 



V x [In f t (x, v ) + 1] • vf t (x, v)dxdv 



(4.68) 



XxB 1 



V v [In f t (x, v) + 1] • [M(x, v, f t )ft{x, v)) dxdv 



XxB 1 



XxB x 



V x f t (x, v) ■ vdxdv - / V v ft(x, v) ■ [M(x, v, f t )] dxdv . 



XxB! 
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The first integral gives zero contribution since j XxBl ft{x,v)dxdv = 1 for all t > 0, i.e. 
ft G x Bi). For the second term notice that V„ • [M(x,v, ft)) = —d, therefore 

/ V v f t (x,v)-[M(x,v,ft)]dxdv = - I ft(x,v)V v -[M(x,v,f t )]dxdv (4.69) 

JXxBx JXxB! 



— d f t (x,v)dxdv = d . 

JXxB 1 

We then obtain (I4.65p . To get (I4.66P we proceed in the same way. We need only to 
modify f!4.69j) as 

V v ft(x,v) ■ [M e (x,vJ t )]dxdv = - [ ft(x,v)V v ■ [M(x,v,ft)]dxdv (4.70) 



= d h t (x)fl(x, v)dxdv . 

JXxBx 

□ 

By the above lemma, if h t (-) ^ for t large enough then 

lim H(f t ) = -oo . 

t— >oo 

From this we can deduce that even starting at time t = from a measure which is 
absolutely continuous with respect to Lebesgue measure in X X B\, having therefore 
finite Boltzmann-Vlasov entropy, at infinity the asymptotic measure is singular with 
respect to the Lebesgue one. 

5. Appendix 

5.1. Proof of Lemma 13.101 We remark that ( 13 . 1 5[) and (" 13 . 1 8[) imply 

N 

^6 M (q(t)) = 0, i = l,..,N. (5.1) 

i=i 

We write ajj(q(t)) defined in (I3.16P as 

C 1 d 

a id (q(t)) = dijiq ) + j ds—a id ((l - s)q° + sq(t)). (5.2) 

We can therefore decompose A(q(t)) as 

A(q(t)) = A(q ) + B(q(t)) , (5.3) 
where B(q(t)) = {k,j(ci(t))h} itj=1 ^ N and 

r 1 d 

bij{q(t)) = J ds—a i , j {{l-s)q° N + sq N {t)). (5.4) 
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N 



Making use of the isomorphism T introduced in (13.231) . it is enough to consider the 
matrix B (q (i)) : R N — > R N since 

||B(q(f))||<||B(q(f)) 

Therefore, setting 



i=l,. .,N ^— ' 



x 



i,3 * 



j{s, t) := (1 - s)(q? - g°) + s ($(*) - q^t)) , i,j = 1, .., iV 



we have 
d 



a itj ((l - s)q u + sq(t)) 



d / ^((i-^)(g?-g?) + ^)-gj(*))) 

rfs I EjLi U r ((1 - s)(q° - ql) + s( qi (t) - q k (t))) 



(5.5) 



(5.6) 



(5.7) 



d 



U R (x itj (s,t)) 



ds VEtiM^M)), 

Vjfrfaifoi)) ■ [-(g° - qj) + g(g - 

£jLitffl(zi,k(M)) 

M^M)) EjLi VC^£x ilfc (s,t)) ■ [-(g? - g fe °) + *(t) - ? fc (t)] 



Efc=i^Kfe(s,t)) 

Hence, since by definition of Ur, sup xgR d < Ur (0) , Vz = 1, .., N, 



N 



Ef=i |V^(x M ( s , 0)| |-(g? - g?) + ?,(*) - q 3 (t) 



J2\hMt))\<2 
< 2 



Ef=i^Kfc(s,t)) 



(5.J 



— ^— sup |W/j(x)| max - g°) + q^t) - qj{t)\ . 



5.2. Proof of Theorem 14.51 We follow [Sj Theorem 5.1 and divide the proof in two 
steps. 

Step 1 We start proving f!4.36j) . Assume that v t and /i t solve (15. 28 p . We have, by the 
triangular inequality, that 

dbcfa, fJ>t) = d bC (vo ° To it [u.],fxo o 7o )t [jU.]) (5.9) 
< d bC (fio o T 0jt [v], fJL o T 0)t \p.]) + d bC (fi,o ° T 0jt [v], v$ o T 0tt [u.]) . 

Denote by w = (x,v),V(fi.) s (w) = (v(s), A(x(s), fi s ) ~ v ( s )) the vector field on the right 
hand side of (14.32)) . The second term can be bounded as 



4c(a*o o T t[u], u o T 0tt [u.]) = e sup 

fev 



TdXB! 



[dH -dp ] (e- Lt foT t , [v.]) 



(5.10) 



<e d fe£ (/io,^o) 
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where L is the Lipschitz constant of V(/x.) s (-). Notice that the Lipschitz bound of 
V(fi.) s (-) can be easily derived from the Lipschitz bound of A(-,/jl.). We get (I5.10P if 
we can show, since f E V, that e~ Lt f o T t) o[iA] is Lipschitz continuous with constant 
one and therefore it belongs to V. Let w(t) = (x(t),v(t)) be the solution of (14. 32ft with 
initial condition wo = (xo, t>o) and let w(t) be the solution of (I4.32p with initial condition 
^o = (xo,vq), then we need to show that 



\f(w(t))-f(w(t))\<C(t)\w -w \ 



with C(t) < e Lt . Writing 



and 



since / E T>, we have 



Furthermore, 



w{t) = w + V(ji.) a (w(s)) 
Jo 



w(t) =w + V(fi.) s (w(s)) 
Jo 



\f(w(t))-f(w(t))\<\w(t)-w(t)\. 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



\w(t) - w(t)\ < \w -w \+ / \VQi.) a (w(s)) - V(tM.) s {w{s))\ds (5.15) 

Jo 

< \wq — wo\ + L I \w(s) — w(s)\ds . 
Jo 

By the Gronwall's inequality 

\w(t)-w(t)\<e Lt \w -w \ (5.16) 

proving e~ Lt f o T 4j o[^.] E T> and so (I5.10p . We are then left with the estimate the other 
term in (15. 9p which, since f ET>, 



4c (/io ° T Q t[v], no o T M) = sup 

fev 



< 



dp {/ o T tfs [v] - f o r t|Q [//.]} 
Ho (dw) \T t! o[v]w - T tfi [n.]w\ =: \(t) 



(5.17) 
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where T t)0 [V.] and T tt0 [n] are both solutions of the equation (I4.32[) with the same initial 
conditions but with different vector fields. We have 



X(t) = / Ho(dw) \{T tfi [v.]w - T t)0 [/i.H} 

JTdXB! 

ft rt 



u (dw) 



T D xB! 

< / Ho{dw) 
+ / uo(dw) 



dsV(v.) s (T sfi [v.]w) - / dsV(ii.) s (T sfi [ii.]w 











ds {V{jy.) s {T Si0 [jy]w) - V(u.) s (T sfi [a.)w)} 



ds{V(ji.) a (Ts,o\pt]w) - V(u.) a (T.Mw)} 



The first term of (15. 18[) can be estimated by the Lipschtiz property of the vector field 



/ fio{dw) / ds{V(v.) s (T s>0 [ 

Jt d xB 1 JO 



v]w) - V(v.) s (T Si0 [n]w)} 



(5.19) 



<L u (dw) ds \T s>0 [u.]w - T afi \p,.]w]\ 

Jt d xb 1 Jo 

= L ds u (dw)\T sfi [v.}w -T sfi [u]w}\ = L / X(s)ds . 

Jo JtdxBi Jo 



For the second term of ( 15. 18ft we have 



fJ>o{dw) 



T D xBx 



ds {V(fi.) a (T a)0 \pi.]w) - V(v.) s (T sfl [n.]w)} 



o 



(5.20) 



< / fJio(dw) / ds\V{fj,.) a {T afi \jj,]w) - V(u.) s (T St0 \p,.]w)\ 
'TdxBx Jo 

ds / fi (dw) \V(fi.) s (T sfi [fi.]w) - V(u.) s (T S fi\p.]w)\ 
>o Jt d xb v 

ds \ n a {dw) \V(/jl.) s (w) - V(u.) 8 {w)\ . 

'0 JTdxBx 



But, 



\V{fi.) s {w)-V{p.) s {w)\ < \A(x,(i a ) -A(x,v a )\ 



(5.21) 



< 



+ 



I Td xb 1 U r( x ~ y) u Vs(dy, du) - J TdxBi U r (x - y)uu s (dy, du) 



! ToxBl U r( x - v)Vs(dy, du) 
It d xb 1 U r( x - V)^(dy, du) - J TdxBi U r (x - y)v s (dy, du) 



I Td xb 1 U r( x ~ V)^s(dy, du) 
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Since for any measure v G M, / rDXBl Ur(x — y)v s (dy, du) > m.i xe j- D Ur(x) = a, we have 
I TdxBi U r( x ~ y)Vs(dy, du) - J TdxBi U r (x - y)v s (dy, du) _ ^ 



< 



I Td xb 1 U r( x - y)^s(dy, du) 



and 



Ir D xB 1 U r( x ~ V) u ^s(dy, du) - J TdxBi U r (x - y)uv s (dy, du) 



d 

i=i 



I TdxBi U r( x ~ V)^(dy, du) 
/tdxBx U r( x ~ y) u i^s{dy, du) - J TdxBi U r (x - y)uiU s (dy, du) 



(5.23) 



< d 



I Td xb 1 U r( x ~ y)Vs(dy, du) 

su PxeT \^ U n( x ) \ + su P*eT D U r( x ) a < \ 
d b c{Us,v s ) ■ 



Therefore, 



\V{u.) s {w)-V{v.) s {w)\<2d 



su Pxer D \VU R (x)\ + sup xeTo U R (x) 



Co 



dbc(Vs,v s ) = —dbc(u s ,K 
a a 

(5.24) 

where we have set c := 2d(swp xe7 - \ VUr(x) | +sup a . gTo Ur(x)). It is essential that a > 0. 
This is the case for interactions considered in the Lemmat a 14 . 1 1 and 14 . 3 1 once the system 
is confined on the torus TiE- Thus, by ( I5U8|) . fl5TT9|) . fl5T20|) and (l5T21j) we have that 



Mt) < L I \(s)ds + — d b c(Hs,v s )ds . 
Jo a Jo 

Hence, since by (I5.18P A (0) = we obtain 



Co 



i L{t - s) d b c(u s ,v s )d& 



X(t)< 

a 

Taking in account (IBTTOj) . (ISTTTj) and ( 15T26|) we get 

dbc{"t, Ht) < e Lt d bC {Ho, ^o) + — / e L(t ~ s] ] d b c{u s , v s )ds . 

a Jo 

Applying the Gronwall's lemma we get bound (14. 36ft . 



(5.25) 



(5.26) 



(5.27) 



Step 2 To prove the existence of a solution for the fixed point equation 

u t = u o T 0jt \ji] , (5.28) 

In the case where Ur is with compact support and M is replaced by M e we have that 

inf (Ur * v)(x) + e > inf U R {x) + e > e. 

In this case X can be either M. d or Td- 
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we use the Banach fixed point theorem. Let \x be the initial condition. To every curve 
t — >■ /i t , /i — /i we obtain the solution curve 

t-^fioToM (5.29) 

Let us denote this map T : Cm — > Cm, where Cm is the space of weakly continuous 
function [0,T] — > M with fi = fi. We equip Cm with the metric 

dM-),u(-)) = sup [e- at d bC (v t ,fi t )] , (5.30) 
te[o,T] 

for some a > which will be suitably chosen. Since (A4, d b c) is a complete metric space, 
so is (Cm, d a )- Now from Step 1 we have 

d b c(vt,lM) = d bC (F(ii(-))(t),r(v(-))(t)) < - f e L «- a) d hC {n»Vs)ds (5.31) 

a J 

cinci therefore 

d a (^(-))(t),T(v(-))(t)) < —^—d a ( A/ (.), !/(.)) (5.32) 

a (a — -LJ 

for a > L. By a suitable choice of a this proves that J 7 is a contraction. □ 
We recall the integral form of the Gronwall's inequality. 

Lemma 5.1. Let L and T positive numbers, f and rj in C([0, T], M). If, for all t e [0, T], 

77(i) < L I V (s)ds + fit) , (5.33) 
Jo 



then 



V(t) < f(t) + L [ e L ^f(s)ds . (5.34) 
Jo 
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